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ABSTRACT 

An increasing amount of recent observational evidence shows that less massive galaxies 
are, the higher on average is their specific star formation rates (SSFR = SFR/Af*, Af* is 
the stellar mass). Such a trend, called the 'SSFR downsizing' (SSFR-DS) phenomenon, is 
seen for local and high-redshift (back to z ^ 1 — 2) galaxy samples. We use observational 
data related only to disc galaxies and explore how does the average SSFR change with z for 
different masses. For all the masses in the range ^ 10^'^ — 10^"'^ M©, the SSFR increases 
with (1 + z) to a power that seems not to depend on M^,, and at all redshifts smaller galaxies 
have ever higher SSFRs; galaxies less massive than Af* ~ 10^° Mq are forming stars at a 
greater rate than in the past assuming constant SFRs over a Hubble time to build stellar mass. 
We show that these features strongly disagree with the A Cold Dark Matter (ACDM) halo 
hierarchical mass accretion rates. Further, by means of self-consistent models of disc galaxy 
evolution inside growing ACDM halos, the effects that disc feedback-driven outflows and gas 
re-accretion have on the galaxy SSFR histories are explored. The parameters of the outflow 
and re-accretion schemes are tuned to reproduce the present-day A/h-A/* relation (A/h is the 
halo mass) inferred from the observationally-based Af* function of disc galaxies. In the case 
of outflows only, the SSFR of individual model galaxies increases with z roughly as (1 + z)^'^ 
for all the masses (somewhat shallower than observations) with a normalization factor that 
depends on mass as A/"^, i.e more massive galaxies have slightly larger SSFRs, contrary to 
the observed strong SSFR-DS trend. For the re-accretion cases, the dependence on z remains 
approximately the same as without gas re-infall, but the correlation on mass increases even 
for most of the reasonable values of the model parameters. The comparison of models and 
observations in the SSFR- A/* plane at z (where the data are more reliable), clearly shows 
the divergent trend in SSFR as the masses are lower (upsizing vs downsizing). We explain 
why our models show the reported trends, and conclude that the SSFR-DS phenomenon for 
low-mass galaxies poses a sharp challenge for ACDM-based disc galaxy evolution models. 

Key words: cosmology: theory — galaxies: evolution — galaxies: haloes — galaxies: high- 
redshift — galaxies: spiral — galaxies: star formation 



1 INTRODUCTION 

The inference of the assembling history of stellar populations in 
galaxies as a function of their morphological type and luminosity 
(mass) is a major topic of present-day extragalactic astronomy, as 
well as a key probe for cosmologically-based models of galaxy 
formation and evolution. 

Based on chemical and spectro-photometric studies of local 
galaxies, it has long been known that more massive/earlier type 
galaxies hold on average older stellar populations and were formed 
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over a shorter time-span, than less massive/later type galaxies (e.g., 
Faber et al. 1992; Worthey et al. 1992; CaroUo, Danziger & Buson 
1993; Bell & de Jong 2000). With the advent of large and com- 
plete surveys of local galaxies and the significant improvement of 
the population synthesis models, the average ages, and even the 
whole star formation rate (SFR) histories of all kinds of galaxies, 
have been inferred (e.g., Kauffman et al. 2003; Heavens et al. 2004; 
Jimenez et al. 2005; Gallazzi et al. 2005,2008; Cid Fernandez et al. 
2005; Panter et al. 2007). The general conclusion of these works 
confirms and extends the previous results: more massive galaxies 
assembled their stars at earlier epochs and on shorter time scales, 
halting their star formation (SF) at later epochs. This phenomenon 
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has been now dubbed as 'archaeological downsizing' (Thomas et 
al. 2005; Fontanot et al. 2009). 

In the last decade, the observational capabilities allowed the 
extension of galaxy population studies to high redshifts. With these 
'look-back studies', the current properties of galaxies at their ob- 
served times, as well as their past histories, were inferred. In a pi- 
oneering work, Cowie et al. (1996) have found that the maximum 
rest-frame K—haad luminosity of actively star-forming galaxies 
has declined with time in the redshift range z = 1.7 — 0.2, i.e. 
the SF efficiency shifts to lower mass galaxies with time. This phe- 
nomenon has been confirmed by subsequent works based on wider 
surveys with multi-wavelenght information (Juneau et al. 2005; 
Feulner et al 2005; Bell et al. 2007; Noeske et al. 2007a,b; Zheng 
et al. 2007; Cowie & Berger 2008; Mobasher et al. 2009; Chen et 
al. 2009; Damen et al. 2009; see Fontanot et al. 2009 for more ref- 
erences). These works allowed to infer the specific SF rates (SSFR 
= SFR/M*, M« is the galaxy stellar mass) for relatively complete 
samples of galaxies at different redshifts until z ^ 1 — 2. The gen- 
eral result is that the measured SSFRs tend to be higher as smaller 
is Af,, a phenomenon called 'SSFR downsizing' , this trend being 
observed at all the studied redshifts. Recent look-back studies of 
galaxy luminosity functions have also confirmed the local infer- 
ences of archaeological downsizing mentioned above: the high- 
mass end of the galaxy mass function seems to be mostly in place 
since z ~ 2, while the abundances of galaxies of smaller stellar 
masses grow gradually with time (Daddi et al. 2004,2007; Bundy et 
al. 2004,2006; Drory et al. 2005; Conselice et al. 2007; Marchesini 
et al. 2009; Perez-Gonzalez et al. 2008 and see more references 
therein). 

The emerging observational picture of downsizing in the stel- 
lar mass galaxy assembly has been confronted with the hierarchical 
clustering scenario of galaxy formation and evolution, based on the 
A Cold Dark Matter (ACDM) cosmological model. At this point, 
as shown in Neistein et al. (2006; see also Fontanot et al. 2009), it 
is important to realize that the downsizing effect to which obser- 
vations refer generically, has actually many manifestations, related 
to different phenomena, involving different types of galaxies and 
different epochs in their histories. Herein we emphasize at least the 
two distinct downsizing phenomena mentioned above: 

(i) the 'archaeological downsizing' related to the early {z^2) 
assembly of most stars in massive/early type, galaxies, and 

(ii) the 'SSFR downsizing' (hereafter SSFR-DS) related to the 
later growth of the relative stellar mass as less massive are the 
galaxies. 

The archaeological downsizing has a partial explanation in the 
frame of the hierarchical clustering process of dark-matter halos 
(Mouri & Taniguchi 2005; Neistein et al. 2006; Guo & White 2008; 
Keres et al. 2009), and astrophysical processes like Active Galac- 
tic Nuclei (AGN)-feedback have been introduced also to get better 
agreement with the observations (e.g.. Bower et al. 2006; Croton et 
al. 2006; Cattaneo et al. 2006; Monaco, Fontanot & Taffoni 2007; 
Hopkins et al. 2008; Somerville et al. 2008, and see more refer- 
ences therein). 

Regarding the SSFR-DS, which is the focus of this work, if 
confirmed, it will turn as a pervasive problem for current mod- 
els. This problem refers mainly to the high SSFRs of low-mass 
(Mt ~ 10^ — 3 X 10^°Mo) star-forming disc galaxies. For these 
galaxies the AGN feedback is not expected to be important because 
either they do not have AGNs (e.g., Kauffmann et al. 2003; Salim et 
al. 2007) or their AGNs are too weak as to produce significant feed- 
back (e.g.. Bower et al. 2006; Crotton et al. 2006). By using semi- 



analytical models (SAMs) several authors have found systemati- 
cally that the SF in modeled low-mass galaxies happens too early 
and is over-quenched at later times, showing these galaxies have 
too low SSFRs at low redshifts as compared with observational in- 
ferences (see for recent results and more references Somerville et 
al. 2008 and Fontanot et al. 2009). The problem is sharpened by 
the fact that in the SAMs, the baryon mass fraction of low-mass 
galaxies has to be decreased systematically in order to reproduce 
the low-luminosity end of the luminosity function (e.g., Benson et 
al. 2003). This is accomplished by introducing appropriately tuned 
schemes of very strong stellar ejective feedback as a function of 
mass, which worsens the SSFR-DS problem (§4). Late re-infall of 
the ejected gas (Bertone, De Lucia & Thomas 2007; Oppenheimer 
& Dave 2008) could offer a partial cure to this problem. 

In this work we aim to 'isolate' the SSFR-DS problem for 
only low-mass disc galaxies and explore in a ttansparent way the 
effects of ejective stellar feedback and gas re-accretion on the 
SSFR histories of galaxies by means of evolutionary hydrodynami- 
cal models of disc galaxies formed inside growing CDM halos. We 
will show that a solution to the SSFR-DS problem is not simple in 
the context of the hierarchical ACDM scenario. 

In Section 2, the observational data to be used in this paper 
are presented and discussed. Preliminary theoretical predictions re- 
garding the SSFR histories of galaxies as a function of mass are 
presented in Section 3. A brief review of our galaxy disc evolu- 
tionary models and the parametric schemes introduced for model- 
ing ejective feedback and gas re-infall are given in Section 4. The 
model results and their comparison to observations are presented 
in Section 5. Finally, a summary of the results and a discussion are 
presented in Section 6. 

The cosmological model used throughout this paper is the con- 
cordance one with h = 0.7, Q.k = 0.7, ^Irn = 0.3, fit = 0.042, 
and as = 0.8. 



2 THE OBSERVATIONS 

Our aim is to target the potential problem of SSFR-DS, and then 
compare observations with detailed theoretical predictions. We fo- 
cus our study here only on normal (not dwarf) low-mass disc galax- 
ies (10^ "' M, < lO^'^"'' Mq). On one hand, due to completeness is- 
sues, dwarf galaxies are not taken into account in the 2; ~ and 
high-redshift observational samples to be used here. On the other 
hand, since most of dwarf galaxies are satellites in bigger systems, 
the physics of them is affected by several enviroimiental processes 
(ram pressure, tidal stripping, interactions, etc.), which are not con- 
sidered in our models. 

In order to determine the SSFR-M, relation at different red- 
shifts, multi-wavclcngth photometric and spectroscopic extensive 
samples are necessary. The stellar masses and SFRs of galaxies are 
inferred commonly by using fitting techniques to stellar population 
synthesis (SPS) models; the SFRs of galaxies can be inferred also 
by using emission lines that ttace SF when available. The uncertain- 
ties associated to these inferences arc large due to selection effects, 
sample incompleteness, and the methods used to infer M« and the 
SFR (see for a discussion §§6.1). For our work, we require samples 
where late-type galaxies are separated from early-t5'pe ones. This 
is the case of the observational works by Salim et al. (2007; « ~ 0) 
and by Bell et al. (2007; high redshifts). 

Bell et al. (2007; see also Zheng et al. 2007 and Damen et al. 
2009) analyze the COMB 0-1 7 photometric redshift survey com- 
bined with Spitzer 24/Lim data (0.2 < zi, 1.0). Stellar masses 
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Figure 1. Evolution of the SSFR of late-type galaxies for different fixed 
values of M». The data at 2 ~ (squared dots) are from Salim et al. 
(2007). The solid curves were infeired from Bell et al. (2007) for their 
two fields, CDFS and A901 (corresponding upper and lower curves, re- 
spectively; see text). From top to bottom, the data correspond to mass bins 
centered at Log(Af*/M0) = 9.5 (cyan), 10.0 (blue), 10.5 (green), and 11.0 
(red), respectively. The thin dot-long-dashed green line is our estimate 
from the data for the individual evolution of an average galaxy that ends 
at z = with Mt = lO^^'^Mo (see text). The solid black line shows 
the curve l/\tu{z) - 1 Gyr](l-R) corresponding to constant SFR along the 
time. The dashed curves are 'dark-matter halo' SSFR histories calculated as 
Mh/[A^h (1 — R)] (see text); from top to bottom the corresponding masses 
areMh = W^-'^Mq (red), = IQI^ IMq (green), A4 = IO^'^Mq 
(blue), and M-^ = 10^^'^Mq (cyan). The corresponding stellar masses are 
roughly as those of the data. The cosmological specific infall rates drive 
SSFRs that ai'e slightly higher for massive galaxies and strongly lower for 
low-massive galaxies than observational inferences. 

were estimated using 17 pass-bands in conjunction with a non- 
evolving template library derived from the Pegase SPS model, and 
SFRs were determined from the combined UV and IR fluxes of 
galaxies. In Fig. [T] we reproduce the averaged data presented in 
their Fig. 3 (SSFR vs. A/,), where four redshift bins centered on 
z ~ 0.9, 0.7, 0.5, and 0.3 are given. We reproduce the data asso- 
ciated only with the blue cloud (late-type) galaxies and for four 
masses [Log(A4",/MQ) = 9.5, 10.0, 10.5, and 11.0 from top to 
bottom, respectively]. Their data are separated according to two 
COMBO-17 fields observed by Spitzer. the extended Chandra Deep 
Field South (CDFS), and the field around the Abell 901/902 (A901) 
galaxy cluster (upper and lower curves in Fig.[T] respectively). Note 
that the SSFRs from CDFS are systematically higher than those 
around the denser field of a cluster. It is well known that the SF 
activity tends to be quenched in galaxies in high-density regions. 
Thus, for comparison with our theoretical predictions, the data cor- 
responding to CDFS are likely more appropriate. 

The data in Fig. [T] corresponding to 2 ~ were taken from 
Salim et al. (2007), who obtained A/* and dust-corrected SFRs for 
^ 50, 000 galaxies by fitting the SDSS and GALEX photometry 
to a library of dust-attenuated SPS models. They were able to sep- 
arate from their sample the "pure" star-forming galaxies with no 



AGN, which form a well defined linear sequence in the SSFR vs. 
A/f, plot, fitted linearly by Log SSFR = -0.35(log M,- 10) - 9.83. 
In Fig. [T] the same four masses related to the data from Bell et al. 
(2007) are used for the Salim et al. (2007) results at 2: ~ 0. 

From Fig.[T]we see that the SSFR of what roughly can be con- 
sidered late-type galaxies declines from 2; ~ 1 to 2 ~ almost 
with the same slope for all the masses. However, the normalization 
strongly depends on A/,, being on average the SSFRs at all red- 
shifts higher as lower is Al,. This shows the SSFR-DS behavior. A 
crude approximation to the range of the data displayed in Fig.[T]is: 

SSFR = 0.13 M-iV"* (1 + zf''^ Gyr~\ (1) 

where Mt,io is M« in units of 10^°Mq. Recently, for all (late- 
and early-type) galaxies and back to 2 « 2.8, Damen et al. (2009) 
have reported a trend of SSFR on redshift also not strongly depend- 
ing on M,, with a trend roughly proportional to (1 + 2)"" (see also 
Feulner et al. 2005; Zheng et al. 2007; Martin et al. 2007; Perez- 
Gonzalez et al. 2008). There are some pieces of evidence that the 
SSFRs of massive galaxies start to increase significantly and over- 
comes the SSFRs of less massive galaxies only at redshifts larger 
than 2 (Perez-Gonzalez et al. 2008) 

Notice that each of the pairs of curves inferred from observa- 
tions for 2> 0.3 in Fig. [T] refer to a given A/, that is the same at 
each redshift, i.e. these curves do not refer to the evolution of in- 
dividual galaxies (evolutionary tracks), but to galaxy populations 
at different epochs. We have carried out an exercise for recover- 
ing approximate galaxy evolutionary tracks from the observational 
data in diagram of Fig.[T] We assume that a galaxy of mass Af* at 
redshift 2 increases its mass during the time interval {dt/dz)dz by 
forming stars at a rate Af, = {1- R)SSFR{M,, z)M,{z), where 
SSFR is given by eq. l[T), and R — 0.4 takes into account for gas 
return due to stellar mass loss. This way, we are able to calculate 
Af,/Af, at each 2 for an "average" galaxy evolutionary track that 
ends at 2 ~ with a given M,. The dot-dashed line shows the 
average evolutionary track for a present-day galaxy stellar mass 
of lO^^'^A^Q. The evolutionary tracks in the range of the data are 
roughly parallel and at 2 ~ intercept the SSFR axis at the level of 
the corresponding masses. As seen in Fig.[T] the evolutionary track 
inferred this way is only slightly steeper than the observed constant 
mass curves. 



3 PRELIMINARY THEORETICAL PREDICTIONS 

From a theoretical point of view, in order to estimate the level of 
SF activity of galaxies as a function of their masses and epochs, 
let consider first the case of a constant SFR history. Then, Af* 
at the cosmic time tH(z) should be given by Af» ~ SFR x(l 
- R)x(t£f(2) - 1 Gyr); I Gyr is subtracted in order to take into 
account the (average) delay in the formation of halos that host 
galaxies. If the measured SSFR at 2 = is smaller (larger) than 
~ l/[tH{z) - 1 Gyr] (1-R), then the average SFR of the given 
galaxy has been higher (lower) in the past than in the current 
epoch. In Fig. [T] the curve l/[t_r/(2) - 1 Gyr] (1-R) for the cosmol- 
ogy used here and for R=0.4 is plotted (thin solid line). The con- 
stant SFR case is just indicative of the situation: low-mass galaxies 
(Af* < lO^^Mo) show current SFRs higher than their past average 
SFR. 

Now, within the context of the ACDM cosmogony, one may 
anticipate a potential disagreement with observations related to the 
SSFR-DS phenomenon. By using an extended Press-Schechter ap- 
proach (Avila-Reese et al. 1998; Firmani & Avila-Reese 2000 - 
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hereafter FA2000; see §4), we calculate for a given present-day 
Mh its averaged mass aggregation history (MAH), My^^z), from 
tens of thousands of Monte-Carlo extractions. By assuming that 
the SSFR is driven by the halo specific mass aggregation rate, 
(Mh/Mh)(2), and correcting by the gas return factor (1 - R), we 
have calculated the 'dark matter'-driven SSFR histories for differ- 
ent halo masses. In Fig.[T]the dashed lines show from top to bottom 
the SSFR histories for halo masses from Mh = IO^^ '^Mq (red) 
to Mh — lO^^'^M© (cyan), which correspond roughly to stellar 
masses from M, = lO^^M© to M, = lO^'^M©, respectively (see 
§§4.1 and Fig. 2). A good interpolating formula to our results since 
2 ~ 1 is: 
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where Af*,i2 is M* in units of lO^^M©. 

From this very preliminary calculation, we conclude that 
while the dark matter-driven SSFR histories of massive galaxies 
do not show a significant discrepancy with observations, for low- 
mass galaxies (M« < lO'^^'^M©) the discrepancy becomes signifi- 
cant and it increases as smaller is the mass. The predicted SSFRs 
of low-mass galaxies are lower than the observed ones beyond the 
uncertainty level. A first guess to solve such a conflict within the 
context of the hierarchical cosmogony is to propose that the main 
stellar formation epoch of low mass galaxies is delayed by some as- 
trophysical processes to lower redshifts. Can gas outflows and later 
re-accretion work in this direction? In what follows, we turn out 
to self-consistent models of disk galaxy formation and evolution to 
explore this question. 



4 THE MODEL 

The formation and evolution of disc galaxies within growing dark 
matter halos can be followed in a transparent and self-consistent 
way by using simplified hydrodynamic models of discs in cen- 
trifugal and vertical hydrostatic equilibrium, and with a SF mecha- 
nism triggered by a disc instability criterion and self-regulated by 
an energy balance process (the semi-numerical approach; FA2000; 
Avila-Reese & Firmani 2000; see for similar approaches van den 
Bosch 2000; Naab & Ostriker 2006; Stringer & Benson 2007; Dut- 
ton & van den Bosch 2009 -hereafter DvdB09). 

The main physical ingredients of the models used here are as 
follows. A special extended Press-Schechter approach based on the 
conditional probability (Lacey & Cole 1993) is used to generate the 
halo mass accretion histories (MAHs) from the primordial Gaus- 
sian density fluctuation field. A generalized secondary infall model 
with elliptical orbits is applied to calculate the time-by-time viri- 
alization of the accreting mass shells (Avila-Reese et al. 1998). The 
orbit ellipticity parameter is fixed in such a way that the structure of 
the ACDM haloes agrees with results from cosmological N-body 
simulations (Avila-Reese et al. 1999; FA2000). A (baryon) fraction 
of the mass of each accreting shell is assumed to cool down in a 
dynamical time and form a disc laye^^. We call 'primary' cosmo- 
logical gas accretion that one related to these accreting mass shells 
with an universal baryon fraction /b.univ = CIb/^dm = 0.163 



^ For halo masses lower than ~ 5 10^^ Mq the radiative cooling time is 
ever shorter than the Hubble time, t[j. Therefore, this assumption works 
well for halos of this mass and smaller. For larger masses, the cooling time 
becomes larger than t[f. Then, the mass baryon fraction available to form 
the galaxy decreases systematically with mass. 



for the cosmology adopted here. The final galaxy baryon fraction, 
/b.gai = Mgai/Mh, is determined by the further disc gas ejecting 
and infalling processes. 

The accreting mass shells at the time of their virialization, tv , 
are assumed to rotate rigidly with a specific angular momentum 
calculated as jsh{tv) = AJ^/AMh, where A represents a differ- 
ence between two time steps, and Jh = XhGM^^'^ / \Eh\^^'^, Mh, 
and Eh, are the halo total angular momentum, mass, and energy 
respectively; is the halo spin parameter, assumed to be con- 
stant in time. As a result of the assembling of these mass shells, 
a present-day halo ends with an angular momentum distribution 
close to the (universal) distribution measured in N-body simula- 
tions (Bullock et al. 2001). The radial mass distribution of the layer 
is calculated by equating its specific angular momentum to that of 
its final circular orbit in centrifugal equilibrium (detailed angular 
momentum conservation). The superposition of these layers form 
a gaseous disc, which tends to be steeper in the centre and flatter 
at the periphery than the exponential law. The gravitational interac- 
tion of disc and inner halo during their assembly is calculated using 
the adiabatic invariance formalism. 

A further step is the calculation of the stellar surface density 
profile. The previously mentioned processes and the fact that SF 
is less efficient at the periphery than in the centre produces stellar 
discs with a nearly exponential surface density distribution, and the 
size mainly determined by Xh- The disc SF at a given radius (as- 
suming azimuthal symmetry) is triggered by the Toomre gas grav- 
itational instability criterion and self-regulated by a balance be- 
tween the energy input due to SNe and the turbulent energy dissipa- 
tion in the ISM. Both ingredients determine the gas disc height and 
the SF rate. This physical prescription naturally yields a Schmidt- 
Kennicutt-like law. The SF efficiency depends on the gas surface 
density determined mainly by A^, and on the gas infall history, 
which in its primary phase is proportional to the halo MAH. 



4.1 Feedback-driven outflows 

In our previous works, we did not take into account disc mass out- 
flows (galactic super-winds) due to SN and radiation pressure feed- 
back. Therefore, the driving factor for disc growth was solely the 
gas infall rate assumed proportional to the halo mass aggregation 
rate, where the proportionality coefficient is the galaxy baryon mass 
fraction, /b.gai. This fraction was fixed to a constant value smaller 
than the universal baryon fraction, /b.univ, by a factor ~ 4. With 
this simple assumption, nearly flat rotation curves at the present and 
at higher redshifts, and good agreement with observed disc scaling 
relations were obtained (e.g., FAOO; Avila-Reese et al. 2008; Fir- 
mani & Avila-Reese 2009). 

Here we introduce an algorithm for outflows from the disc 
similar to those used in the SAMs but applied locally to the evolv- 
ing disc-halo system as a function of the radius (DvdB09). The 
outflow is assumed to move at the local escape velocity of the disc- 
halo system, Vcsc(''). Therefore, the efficiency of mass ejection de- 
pends on mass. For halos less massive than ~ 10^^ Mq, the out- 
flows were found to be able to reduce the initial universal baryon 
fraction by factors ~ 2 and more as the mass is smaller (c.f. van 
den Bosch 2002; Benson et al. 2003; DvdB09). The mass ejected 
per unit of area and time from disc radius r is given by 



t(r) = EsFR(r-) e 



/ 1000 kms- 
l V;ec(r) 



(3) 



where Esfr is the SFR surface density, and e is a free param- 
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Figure 2. Panel (a): Semi-empirical and modeled stellar vs halo masses. 
The former are inferred by different authors by matching a given observed 
luminosity (Mj,) function with the theoretical halo mass function for all 
galaxies and halos (non-continuous curves), for only late-type galaxies and 
halos that likely host late-type galaxies (continuous red curve for the av- 
erage, and dotted curves encompassing the dotted horizontal lines for the 
1 (T uncertainty), and from direct weak gravitational lensing studies (cyan 
solid squares with vertical error bars). The con'esponding literature sources 
are indicated inside the figure. For the model predictions (thick dot-dashed 
curves), the disc outflow cases for n = 1 , 2, and 5 (respectively in clock- 
wise sense) are shown. The models are forced to agree with the continuous 
red line at 10^" by tuning the free parameter e in eq. (3ji see for the 
values Table 1. The n = 2 'energy-driven' outflow case best reproduces the 
observational inferences. Panel (h): Same as panel (a) but for models that 
include gas re-accretion. Only cases with n = 1 and Afpcg=const. {fi = 0, 
see eq.|5) are plotted. For larger values of n, the re-accretion produces too 
shallow slopes in the M-^^-Mf relation. The three cases plotted with thick 
solid, dashed, and dot-dashed lines, correspond to values for the (Log(m), 
a) parameters of (10.0, 1.0), (10.5, 1.0), and (10.0, 0.5), respectively (see 
Table 1). 

eter. In the literature, values for n from 1 to 5 were commonly 
used in order to reproduce the low-mass luminosity function. The 
values n = 2 and 1 are expected for the 'energy-driven' and 
'momentum-driven' outflows, respectively. For the former case, 
assuming that each SN produces an energy Esn ~ 10^^ erg, 
and that the the number of SNe per solar mass of stars formed 



is 77SN = O.OO9M0~^, one obtains namely the value of ~ 1000 
kms~^ reported in eq. lO and the parameter e can be interpreted 
as the fraction of SN kinetic energy transferred into the outflow. 
For the latter case (n = 1), by assuming that each SN produces 
a momentum psN ~ 3 10* M0kms~^ and using the same value 
of rysN as above, the normalization velocity in eq. ((S} should be 
~ 300 kms~^; by keeping 1000 kms~^ in eq. lO, the outflow mo- 
mentum is already 3.3 times the momentum produced by one SN; 
this should be still multiplied by our e. However, we should say that 
the 'momentum-driven' case is thought for the effect of massive- 
stars rather than for the SN contributions (e.g. Murray, Quataert & 
Thompson 2005; Oppenheimer & Dave 2008; Grimes et al. 2009); 
the total momentum (energy) injected by massive stars during their 
life is much higher than the one injected by SNe. 

Because the outflow rate is proportional to the SFR (eq.O, the 
outflows are more efficient in removing gas at earlier epochs, when 
the SFR is higher. The early accreted gas has lower angular mo- 
mentum than the later accreted gas. Thus, as the result of disc gas 
ejection by outflows, the spin parameter A of the final disc becomes 
larger than the spin in case of no outflows, which, by assumption, 
is equal to the halo spin parameter (see also DvdB09). 

The study of the galaxy-size outflow physics is a complicated 
hydrodynamic and radiation transfer problem, and it requires a very 
large dynamical range, from parsec to megaparsec scales (e.g., 
Dalla Vecchia & Schaye 2008; Oppenheimer & Dave 2008 and 
more references therein). The model described above is certainly 
an oversimplification, which overestimates the feedback efficiency 
(see for a discussion DvdB09). In most of the previous works it 
was assumed that the ejected gas is lost forever by the galaxy- 
halo system and the proportionality coefficient and exponent in eq. 
([3) were varied in order to reproduce the low-luminosity side of 
the observed luminosity function as well as galaxy scaling laws 
like the TuUy-Fisher relation (c.f. Benson et al. 2003; DvdB09). 
Feedback-driven outflows seem to be also necessary to explain the 
steep mass-metallicity relation of galaxies at redshifts z ~ 1 — 2, 
and the observed metallicities of the IGM at high redshifts (Aguirre 
et al. 2001; Oppenheimer & Dave 2006; Finlator & Dave 2008; 
DvdB09). 

The full modeling of feedback-driven outflows is out of the 
scope of this paper. Instead, we aim to explore in a semi-empirical 
way the effects of different kinds and levels of feedback on the 
SFR history of disc galaxies as a function of mass. The exploration 
is semi-empirical in the sense that we take care to agree with the 
present-day J\fh-A/, relation inferred from matching the galaxy 
stellar (luminosity) and halo mass cumulative functions. In Fig. |2] 
we show this relation as inferred by different authors (Shankar et 
al. 2006; Baldry, Glazebrook & Driver 2008; Conroy & Wechsler 
2009; Moster et al. 2009; Rodn'guez-Puebla et al. 2009, in prep.). 
Different observed galaxy luminosity functions, different methods 
to pass from luminosity to A/,, and different halo mass functions 
(analytical fits or halos directly from N-body numerical simula- 
tions), as well as different corrections to account for halo groups 
and sub-halos, were used by each one of these authors. For com- 
pleteness, we reproduce in Fig.|2]also the more direct (but yet very 
uncertain) inferences based on weak lensing studies for late-type 
galaxies (squares with large vertical error bars; Mandelbaum et al. 
2006). 

The M\i-Mf relation has been inferred commonly using the 
whole galaxy population and all the halos. In the case of Baldry 
et al. (2008; long dashed blue line), the used galaxy sample refers 
to galaxies in the field, where disc galaxies dominate. Since our 
study refers to disc galaxies, the A/h~Af, relation should be in- 
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ferred from a stellar mass function for only disc galaxies and from 
a halo mass function related to halos that will host disc galax- 
ies. This inference has been performed in Rodriguez-Puebla et al. 
(2009), who used the A/, function for SDSS late-type (blue) cen- 
tral galaxies as reported in Yang, Mo & van den Bosch (2009), 
and the halo Afh function corrected for halos that did not suffer a 
major merger (mass ratio larger than 0.2) since z = 0.85; Gover- 
nato et al. (2009) have shown that after a major merger took place 
before z ~ 0.8, the galaxy can regenerate a significant disc until 
the present epoch. Rodrfguez-Puebla et al. checked that the frac- 
tions of late-type galaxies and 'quiet' halos with respect to their 
corresponding distribution functions were similar, around 55% in 
both cases. The solid red line in Fig. [2] reproduces the results by 
Rodrfguez-Puebla et al. including the 1 o uncertainty (dotted red 
lines which encompass the area of dotted horizontal lines); this un- 
certainty is mainly due to the uncertainties in the population syn- 
thesis models used to estimate Af, . 

We use the A/h-Af, relation by Puebla-Rodriguez as indica- 
tive, intending that the feedback-driven mass ejection shall help 
reproduce roughly such a (semi-empirically derived) relation. It 
should be stressed that our conclusions below regarding SSFRs are 
unchanged if our models are calibrated to reproduce any of the 
other Mh-M, relations showed in Fig.|2] We focus our study only 
on galaxies with Mt < 10^°-^Mq, those for which the SSFR-DS 
phenomenon becomes strong (see §2). 

In panel (a) of Fig.|2] the thick dot-dashed lines show our re- 
sults for the cases n — 1,2, and 5 in clockwise sense, respectively. 
In each model the outflow efficiency e is fixed in order to obtain 
the best agreement with the A/h-Af* relation at Af* = 10^°Mq. 
The higher n is, the shallower the Afh-A^» relation at low masses; 
for example, the slopes in the range 9.5 ^ LogAf, ^ 10.5 are 0.79, 
0.62, and 0.36 for n = 1,2, and 5, respectively (the correspond- 
ing slope for the Rodrfguez-Puebla et al. inference is 0.65). We 
have experimented also with an outflow model with constant mass 
loading (i.e., n — 0); as expected, the models follow a correla- 
tion steeper than the case n = 1. If outflows are not introduced 
(/b,gai=0.04=const. is used), then the correlation is only slightly 
shallower than the case Afh oc Af, in Fig.|2](slope ~ 0.94); this is 
because in our models the lower the disc mass is, the less efficient 
is the process of gas transformation into stars. 

As seen from Fig.|2] a reasonable agreement below lO^^'^M© 
with the observational inferences is obtained for n = 2 ('SN 
energy-driven' outflows) with a high SN energy transference ef- 
ficiency, e — 0.62. The values of e for different cases, as well as 
the obtained values of A//* at z = for different halo masses, are 
presented in Table 1 . 

4.2 Re-accretion of the ejected gas 

The assumption that gas is lost forever is strong and likely unre- 
alistic. The outflow propagates along the intra-halo and/or intra- 
filament media and even if gas escapes from the current gravita- 
tional potential, it will be slowed down and stopped by further in- 
teraction with the intergalactic medium, and then re-accreted later 
by the halo-galaxy system (see e.g., Bertone et al. 2007; Oppen- 
heimer & Dave 2008). 

The motion of an ejected disc gas element can be calculated by 
taking into account the halo-galaxy gravitational potential and the 
viscosity of the surrounding gas. Once the ejected gas element stops 
with respect to the gas around, both the ejected gas element and the 
circungalactic gas will infall together onto the disc galaxy. This is 
a rather complex hydrodynamic process, demanding in computing 



time and, in any case, uncertain due to the stochastic nature of the 
problem. Instead, we adopt here a simple approach that can be eas- 
ily implemented in our numerical code. Let us consider that during 
the time interval elapsed between the ejection of a gas element and 
its re-infall onto the disc, a certain amount Mpcg of primary cos- 
mological gas (see above for a definition) accretes onto the disc. 
The information related to the surrounding gas that breaks the out- 
flow (mostly related to the infalling cosmological gas), as well as of 
the dynamics in the gravitational field, are actually summarized in 
the parameter A/pcg (see for a similar statement Dubois & Teyssier 
2008). In other words, the physics concerning to the outflow/re- 
accretion process can be parametrized through A/pcg. With this 
idea in mind, we follow the next simple approach for calculating 
gas re-accretion in our code. For a given disc gas element, as soon 
as is ejected, the acumulated amount of the accreted primary cos- 
mological gas is numerically calculated time by time. When such 
amount reaches the value of A/pcg, then the given mass element is 
reintegrated to the disc. This procedure is applied to each ejected 
gas element at any time. Adopting such an approach, the physics 
reduces to fix A/pcg in a parametric fashion. 

Actually we introduce a lognormal distribution around our pa- 
rameter A/pcg in order to take into account the stochastical nature 
of the problem. Such distribution is defined by 

dP = e-"'('"(^^P-g)-<'"(*^P-g)0'd;n (A/^cg) , (4) 

where a regulates the width of the distribution around the loga- 
rithmic mean value (/n(A/pcg)); a and A/pcg= e*<'"'^^pcg'>) are 
introduced as free parameters. 

By means of cosmological hydrodynamical simulations, Op- 
penheimer & Dave (2008) have found that the re-accretion time 

— 1/2 

scales with galaxy mass roughly as A/^^j and they interpret this 
to be the case when environmental effects are dominating the retar- 
dation of outflows. Such a behavior implies in our case that A/pcg 
should be roughly constant or moderately decreasing with maso 
This way, for massive galaxies, A/pcg is attained soon and gas 
re-accretion becomes consequently efficient (galactic fountains), 
while for lower mass galaxies, attaining A/pcg takes longer times; 
for the smallest galaxies, it may takes eventually periods longer 
than the current Hubble time in which case the gas does not re- 
turn. Given the large uncertainties in the physics of the process, we 
introduce a parametric expression for A/pcg: 

where m and fj, are parameters to be probed. 

^ The properties of A/pcg can be estimated also by simple momentum or 
energy conservation arguments. Let us consider a given outflow of mass 
A/o moving in a given solid angle at the ejected velocity Vo in the halo 
potential; then, MoV^ = {Mo + A/p) V°-, where Mp is flie mass piled 
by the outflow, V is the shock velocity, and a = 1 or 2 for momentum or 
energy conservation, respectively. A rough and general estimate for V could 
be the typical halo circular velocity, Vc, which approximately scales with 
1/3 

the halo mass as Vc oc . Therefore, the mass piled by the outflow 

is A/p S3 {MoVo)°-/M'^^^. Such mass will drag the mass element A/o 
back to the disc. In other words, Mo is reintegrated to the disc when, within 
the given solid angle, the amount of mass A^p of primary gas is accreted. 

Thus, A/pcg is just Mp scaled to the entire soUd angle 47r. From this crude 

1/3 2/3 

approximation, A/pcg oc M^ ' or M^ ' , dependencies that are within 
the range explored below by us. 
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Table 1. Predicted A/» and SSFR at 2 = as a function of halo mass Afj^ for different models with only outflows (lines with columns 3-5 empty) and with 
both outflows and gas re-accretion. All the masses are in units of Mq, and the SSFR is in units of Gyr~^. 



Regarding the angular momentum of the re-accreted gas, the 
situation is even more uncertain. We assume that the specific angu- 
lar momentum of this gas is a fraction of the one corresponding 
to the currently hierarchically infalling cosmic gas. The results to 
be presented here depend poorly on the assumed fraction within the 
reasonable range of fj « 0.3 — 0.7 

In panel (b) of Fig. |2l the thick black curves show our re- 
sults for the outflow model with n — 1 (see eq. [3j and constant 
values of the re-accretion Mpcg parameter, i.e. with fi — (see 
eq. |5}: m = 10^° (solid line) and m = lO"^'^ (dashed line), in 
both cases using a = 1 in the probability distribution of i\fpcg 
(eq.|4ll; and m — 10^" but a — 0.5 (dot-dashed line). The corre- 
sponding slopes of the Afh-M, relation in the range 9.5 5; LogA/, 
< 10.5 are 0.58, 0.57, and 0.52, close to the slope of the observa- 
tions. For outflow models with n larger than 1, after taking into 
account re-accretion, the A/h-Af* relations become significantly 
shallower than the observational inferences. In Table 1 we report 
our results at 2: = (Af* and SSFR for different halo masses) for 
the re-accretion parameters mentioned above as well as for some 
cases with /i 7^ 0, and for n = 2. As expected, by including re- 
accretion and trying to reproduce the A/h-A/, relation implies very 
high outflow efficiencies. We see also that the results depend actu- 
ally weakly on the dispersion parameter a (e.g., compare rows 4 
and 6 from Table 1). 



4.3 Strategy 

A given galaxy model is defined by the halo mass, Afh, its MAH, 
and its A^. The initial baryon fraction is assumed equal to the uni- 
versal one, /b,gai=/b,Univ, but the feedback-driven outflow and gas 
re-accretion parameters define the actual value of /b.gai- Since we 
are interested here in generic evolutionary trends related to the SFR 
and M, histories as a function of mass, we study only the "central" 
models of different masses characterized by: 

• the averaged MAH corresponding to the given halo mass Afh, 

• a value of — 0.03, which is somewhat smaller than the 
mean of all relaxed haloes measured in numerical simulations (e.g.. 
Belt et al. 2007). 

Note that these physical ingredients have in fact probability den- 
sity distributions (not taken into account here) that, of course, will 
produce an intrinsic scatter in the values of the galaxy properties 
studied. 



We first study models with galactic outflows only, and then 
include the possibility of gas re-accretion. In all the cases the pa- 
rameters of the feedback-driven outflows and gas re-accretion are 
chosen appropriately to reproduce the A/h-Af, relation inferred 
semi-empirically (Fig. 2). Our aim is to explore whether the model 
SSFR histories for different masses agree or not in general with the 
observational trend of SSFR-DS discussed in §2 (Fig.[Tll. 



5 RESULTS 

The observational data in the SSFR vs. (1 + z) diagram (Fig.[Tll re- 
fer to curves of constant corresponding to galaxy populations 
at different redshifts rather than to curves equivalent to the evolu- 
tion of individual galaxies. However, as discussed in §2, at least for 
z^l, the later curves seem to be only slightly steeper than the for- 
mer ones in the SSFR vs. (1 + z) diagram. Having in mind this small 
difference (note that at z ~ 0, the comparison of models and ob- 
servations is completely fair), we compare below the observational 
data with results for individual galaxy evolution models, which end 
at = with stellar masses similar to those corresponding to the 
data: Log (Af./Mo) = 9.5, 10.0, and 10.5. The massive case. Log 
(Af*/M0) = 1 1.0, as has been explained in §4, is out of the scope of 
the problem studied herein. Recall that our models are constrained 
here to agree with the Afh-Af, relation inferred semi-empirically 
(Fig.H. 

5.1 Galactic outflows only 

We first experiment with the usual ejecting SN feedback model, 
where the gas in the outflows (galactic super-winds) is assumed 
to be lost forever from the halo (§§4.1). In Fig. |2] we showed 
that the semi-empirically derived M^^-M^ relation is better repro- 
duced by the n = 2 outflow model ('energy-driven' outflows). 
The SSFR evolution for this case and for the present-day masses 
Log(A/,/M0) = 9.5 and 10.5 are plotted in Fig. |3] with the dotted 
cyan (lower curve) and green (upper curve) lines, respectively. The 
SSFR tracks corresponding to masses in between this mass range 
lie within these two curves. 

For all the masses, the SSFR increases with z with a trend in 
rough agreement with observations, but there is a small but system- 
atical trend of higher SSFRs as the mass is larger, opposite to the 
strong obsened SSFR-DS phenomenon. A good fit to the model 
results from 2 = to 2 ~ 3 is: 
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Figure 3. Same as in Fig. [T]but including evolutionary tracks of individual 
galaxies. The dotted curves are for models with outflows and for the favorite 
case n = 2; the upper green and lower cyan lines are for Mt = 10^" 
M0 and M, = 10^'^ Mq, respectively. The dashed curves are for models 
with outflows (n = 1) and gas re-accretion = 0, m = 10^", a = 1.0); 
the upper green and lower cyan lines are for = 10^"-^ Mq and M* = 
lO^'S M0, respectively. 



SSFR = 0.005 M,% (1 + zf-"^ Gy^ . (6) 

Such a behavior is consequence of (i) mainly the CDM halo MAHs 
which drive the gas infall process and, consequently, the SF histo- 
ries: the less massive the halos, the lower the specific mass aggre- 
gation rates at late epochs (see Fig. [T); and (ii), in a less extent, of 
the fact that the disc outflow gas ejection is more efficient for lower 
mass halos. The SSFR histories of low-mass models have signif- 
icantly over-quenched present-day SFRs, i.e. their stellar masses 
had to be assembled with SFRs in the past higher than the current 
one. In general, all the models since z ~ 2 lie below the curve 
l/[tH(z) - 1 Gyr] (1-R) of approximately constant SFR history 
(thin black solid line in Fig.[T] see §3). 

The most reliable observations are for local galaxy samples. 
In Fig.|4]we plot the fit and its 1 o dispersion to the SSFR-Af, cor- 
relation given for local SDSS galaxies by Salim et al. (2007). They 
report the fit for only star forming galaxies (solid blue line) and for 
both this sample and the one of galaxies with AGN (short-dashed 
red line). In both cases the galaxy SSFRs significantly decrease 
as M« increases. The three long-dashed magenta curves shown in 
this plot correspond to models with only outflow and with values 
of n = 5, 2, and 1 from top to bottom, respectively (see also Table 
1). Thick line is used for our preferred n — 2 case (see panel (a) of 
Fig.|2]l. The difference between models (weak upsizing) and obser- 
vations (strong downsizing) in the trend of the SSFR-A/, relation 
is significant. In particular, for M» < 1O^°'^M0 , the lower the mass, 
the bigger the differences between the model SSFRs and those in- 
ferred from observations. For A/, > IO^^'^Mq, the differences 
tend to disappear. 



5.2 Galactic outflows + re-accretion 

In §§4.2 we described our parametric model to account for the re- 
accretion of the ejected gas by feedback-driven galactic outflows. 
The models with re-accretion agree with the Mh-A'h relation only 
for outflows with n — 1 (see panel b of Fig.|2j. In Fig.|3] the case 
with n — 1, fi = 0, m — 10^" (see eq. [5) and a — 1 (see eq|4ll 
is plotted for two models that end at z = with Log(Af,/M0) = 
9.5 (upper dashed cyan line) and 10.5 (lower dashed green line). 
We have experimented also with many other re-accretion parame- 
ters (see Table 1). As expected, the SSFRs are higher for the mod- 
els with re-accretion. This is because the rate of later re-accretion 
of the ejected gas sums up to the hierarchical halo accretion rate, 
raising the SFR. From Fig. |3]we see that the SSFR increases with 
(1 + .z) to a power not strongly dependent on mass and actually, 
similar to the case of no re-accretion (dotted curves; see eq. ll6l). 
It is obvious that if re-accretion is included, then the mass load- 
ing in the outflow model should be increased in order to reproduce 
the same Afh-A/* relation than without re-accretion. For example, 
compare the models given in rows 1 and 4 in Table 1 (n = 1 out- 
flow case), where e increased from 1.83 to 7.23; note that the SSFR 
increased for all the masses, but more for the more massive models. 

The physics of the model is rather complex. For A/pcg= 
const., the larger the mass, the shorter the time period of gas re- 
incorporation into the disc because Mpcg is equaled quickly by the 
accreted primary cosmological mass in more massive galaxies. Fur- 
thermore, a given mass element can be ejected and re-accreted sev- 
eral times depending on the outflow efficiency. Given that massive 
galaxies recover their ejected gas sooner than low mass galaxies, 
then the same stellar mass (Afh-Af* relation) is reached decreas- 
ing the power n of eq. lO, i.e. decreasing the low mass galaxies 
outflow with respect to massive galaxies. 

The behavior of the SSFR with mass is again against the 
SSFR-DS. For the value of A/pcg= Mq used here, models 
with A/. ~ 10^° '^ Mq raise their SSFR to values slightly lower 
on average than observational inferences, while the lowest mass 
models, only very slightly raise their SSFR at late epochs, far away 
from the values required to reproduce the high SSFRs inferred for 
low-mass disc galaxies. 

Further, we have experimented with cases where /i in eq. ^ 
is varied. For ^ < and for a fixed value of m, the epoch when 
the SSFR becomes higher due to gas re-accretion is delayed with 
respect to the case /x = 0. This increase is early and large for more 
massive galaxies, but with time the SFR-ejection-re-accretion pro- 
cess is stabilized. For small galaxies, the influence of re-accretion 
over the SSFR begins at late epochs. Overall, taking as an exam- 
ple the case /i = — 1, galaxies that end with Af* > 10^°'^ Mq 
have SSFR histories close to those inferred from observations. For 
smaller galaxies, while their SSFRs are now higher than in previ- 
ous cases, they remain being far from the observational constraints. 
Models keep showing a trend in disagreement with the SSFR-DS 
trend. Regarding the cases with ^ > 0, the SSFRs of less massive 
galaxies result even lower than in the jj, — case. 

In Fig. |4] the loci of the different re-accretion models in the 
z = SSFR- A/, diagram are plotted (see also Table 1). The /i — 
0, m = 10^" fiducial case is shown with the thick dot-dashed black 
line. The uppermost dot-dashed curve corresponds to our /i = — 1 
case. The dot-dashed curve immediately lower to the fiducial case 
is for the same values of this case but for q = 0.5. The next lower 
curve is as the fiducial case but for m = 10^°'^. Then follow the 
curves corresponding to our n — 1 and n — 2 cases (see Table 1). 
As seen, for a reasonable range of values for the re-accretion model 
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Figure 4. Local (z ~ 0) SSFR-Af, relation. The fit to the SDDS data 
(average and 1 a scatter) by Salim et al. (2007) is plotted with solid blue 
lines for only star forming galaxies, and with short dashed red fines for botfi 
the star forming and AGN galaxies. The long-dashed magenta curves are 
for models with only outflows and n = 5, 2 and 1 (respectively from top to 
bottom). Our preferred case, n = 2, is highlighted with bold. The dot-long- 
dashed black curves are for models with outflows and gas re-accretion. Our 
preferred case, n = 2, p. = 0, m = 10^", a = 1.0 is highlighted with 
bold. For the other curves, see text and Table 1. In all the cases, models 
show a weak 'upsizing' behavior contrary to the strong 'downsizing' of the 
observations. 

parameters, the SSFR as a function of Af, remains in conflict with 
observations, and the SSFRs of low-mass models are too low. 



6 SUMMARY AND DISCUSSION 

By means of self-consistent models of disc galaxy evolution inside 
growing ACDM halos, the effects of feedback-driven disc ejective 
outflows (galactic super-winds) and gas re-accretion on the evolu- 
tion of the SSFR for low-mass galaxies (Af* < lO^^'^M© at 2: = 0) 
have been explored. We have studied only the 'central' (average) 
models corresponding to different halo masses since our main goal 
was to probe general trends regarding the evolution of the SSFR 
and its dependence on mass. 

6.1 The effect of gas outflows 

The rate of mass ejection was assumed proportional to the local 
SFR and inversely proportional to the local escape velocity Ksc to 
a given power n (see eq. [3] for n = 2 and n = 1 the outflows are 
expected to be energy- and momentum-driven, respectively). The 
parameters of the outflow models were explored in the light of the 
Mh-M, relation at z = (galaxy efficiency) inferred from a semi- 
empirical approach that makes use of the observed luminosity (A^*) 
function (Fig.[2ll. 

Our best outflow models were chosen to be those that agree 
with the Afh-Af, relation inferred for disc (blue) galaxies and halos 



that did not suffer a major merger since z — 0.85. We have found 
that the slope of the Afh-Af. relation at masses Af* < 3 10^° Af, is 
better attained for the outflow model with n — 2 (energy-driven 
SN feedback, eq.[3} and for a high SN energy efficiency, e — 0.62 
(Fig.m panel a). 

For the n = 2 outflow case, the SSFR of all masses raises 
with redshift proportional to (1 + zf''^ up to z ~ 3 with a little 
dependence on (present-day) mass, as A//."-^ (eqM- The SSFRs 
at low redshifts for models with Af, < IO^'^'^Mq are well below 
the constant SFR curve l/[tH(z) - 1 Gyr](l-R) (see §3), which 
means quiescent late SF activity, and miserably fail in reproducing 
the observational inferences (Fig. [3}. In particular, the models show 
a trend contrary to the observed SSFR-DS phenomenon. Such a 
conflict is clearly seen in the SSFR vs A/, diagram for local galax- 
ies (z ~ 0), where the observational data are more reliable: while 
the SSFR of observed galaxies significantly decreases with A/*, the 
models show a slight increase with A/* . The situation is similar for 
model outflows with n = 1 and n = 5 (Fig.Q. 

We stress that the conflict is not related to the observed low 
SSFRs of massive disc galaxies ( Z IO^^'^Mq) with respect to mod- 
els, but instead is related to (i) the observed trend of SSFR increas- 
ing as the mass decreases (SSFR-DS), and (ii) to the too high val- 
ues of the SSFRs of galaxies with masses M, < 10^'^ Mq, both 
items applying since z ~ 1. 

In the models, on one hand, the SFR history is largely driven 
by the hierarchical halo mass aggregation ratfl which on average 
increases with redshift at nearly the same rate for all masses, but at 
a given z, massive halos have rates slightly higher than those of less 
massive halos (see Fig.[T]and eq.O. On the other hand, the ejective 
feedback scheme used here produces proportionally more gas loss 
(hence, less later SF) for the lower mass discs than for the massive 
ones. We conclude that the low SSFRs of low-mass disc galaxy 
models with outflows, and the (slight) increase of such SSFRs with 
A/*, opposed to the observed strong SSFR-DS phenomenon -the 
lower the mass, the higher the SSFR, is a natural consequence of 
the ACDM halo assembling and, in a minor level, of the feedback- 
driven outflow schemes commonly used to reproduce the local Afh- 
relation. 



6.2 The contribution of gas re-accretion 

We further explored the possibility that the gas ejected from the 
disc-halo system can be re-accreted later onto the disc. Our scheme 
for such a process is very general and encompasses a large range of 
possibilities for the unknown and complex outflow hydrodynamic 
and radiation transfer processes. Instead of introducing as a key 
parameter the time elapsed between when a given mass element 
is ejected and it is further re-accreted, we use as a parameter the 
primary cosmological gas mass that will be accreted during this pe- 
riod, A^pcg- This way, we take into account in some way the amount 
of circumgalactic mass that interacts with the outflow, amount of 
mass that in our models is related to the MAH of the given halo. 
We use a parametrization for A/pcg given by eq. J5](> ™d ^ lognor- 
mal distribution (eq.|4l( for such a parameter is considered in order 
to introduce a (natural) dispersion around the re-accretion times of 
the ejected mass shell. 

In our models, the SFR history depends also on the local disc gas sur- 
face density and on some disc dynamical and hydrodynamical properties. 
Nevertheless, the dominant factor in the SFR histories of our modeled disc 
galaxies is the cosmological gas infall rate. 
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Our results show that the gas re-accretion may raise the SSFR 
of disc galaxies but it does it in the incorrect direction regarding 
mass: while for low mass galaxies the increase in SSFR is small, 
for the massive galaxies the increase is significant at all epochs 
(Fig- Ell. We have experimented with several cases: Mpcg constant 
for all masses, and Mpcg decreasing or increasing with mass. The 
decrease of SSFR as M, decreases (upsizing) is seen in all the cases 
(for 2 ~ 0, see Fig.|4j. 

We conclude that for models of disc galaxy formation and 
evolution in the context of the ACDM cosmogony, the problem of 
the SSFR-DS is not solved by an interplay of outflows and re- 
accretion of gas. 

6.3 Is the SSFR-DS phenomenon well established? 

We have shown in a clear and transparent way the difficulty galaxy 
evolution models in the context of the ACDM cosmology have for 
explaining the SSFR-DS phenomenon, related mainly to sub-L* , 
star-forming (late-type) galaxies. Before discussing some caveats 
of the models and possible cures to the problem, it should be 
stressed that at the level of obsen'ational inferences there are still 
large uncertainties. 

The use of SPS models to fit observational data and infer A/* 
and SFR should be taken still with caution due mainly to the uncer- 
tainties in stellar evolution (for example in the thermally-pulsating 
asymptotic giant branch, TP-AGB, and horizontal branch phases) 
and to our poor knowledge of the initial mass function, IMF, as well 
as due to degeneracies like the one between age and metallicity 
(see for recent extensive discussions Maraston et al. 2006; Bruzual 
2007; Tonini et al. 2009; Conroy, Gunn & White 2009a; Conroy, 
White & Gunn 2009b). For example, Conroy et al. (2009a) esti- 
mate that stellar evolution uncertainties can introduce up to ~ 0.3 
dex statistical error in M» at z ~ 0. Regarding systematical errors, 
they are more difficult to evaluate. Maraston et al. (2006) claim that 
the stellar masses of galaxies with dominating stellar populations 
of ~ 1 Gyr age could be on average 60% lower if their assump- 
tions for convective overshooting during the TP-AGB phases are 
used. It is not easy to estimate the direction in which these statis- 
tical and systematical uncertainties could influence the SSFR-A/, 
dependencies at different redshifts; at least, it is not obvious that 
the SSFR-DS phenomenon could be eliminated. 

It is healthy to stress that the approach we follow for com- 
paring models and observations is the recommended one in the 
sense of minimizing uncertainties. Conroy et al. (2009b) show that 
it is better to first infer the physical properties of observed galaxies 
(e.g., stellar mass and SFR) by using the SPS technique, and then 
compare them to galaxy evolution models, rather than applying the 
SPS technique to models in order to compare their predictions with 
the direct observables. 

Other significant sources of uncertainty in the inferred SSFR- 
M* relations at different redshifts are the selection effects due to 
the incompleteness of the sample, dust absorption, environmental 
effects, and/or limit detections of the tracers of SF due to flux- 
limits or low emission-line signal-to-noise ratios; in addition, ob- 
scured AGN emission could be contaminating the infrared flux and 
some of the optical lines used to estimate SFR (e.g., Daddi et al. 
2007; Chen et al. 2009). 

The main concern regarding the SSFR-DS phenomenon 
among these issues is that selection and environmental effects could 
bias the observed trend that the SSFR increases as A/* is smaller. 
For example, the high SSFRs of low-mass galaxies could be due to 
transient star bursts if the SF regime of these galaxies is dominated 



by episodic processes (external like mergers or internal like statisti- 
cal fluctuations in massive star formation); if among the low-mass 
galaxies, those with low-SFRs are missed due to detection limits, 
then the SSFRs of low-mass galaxies will be biased on average to- 
ward higher values of SSFR, a bias that increases for samples at 
higher redshifts. Nevertheless, in most of the observational stud- 
ies reporting the SSFR-DS phenomenon the authors discuss that, 
while this is possible at some level, hardly would it be the dominant 
factor in the SSFR-DS phenomenon observed since 2 ~ 1 (see for 
example Noeske et al. 2007a). 

On the other hand, the analysis of very local surveys certainly 
helps to disentangle whether episodic star bursting events dominate 
or not the SF history of low-mass disc galaxies. From a study of 
the SF activity of galaxies within the 1 1 Mpc Local Volume, Lee 
et al. (2007; see also Bothwell, Kermicutt & Lee 2009) have found 
that intermediate-luminosity disc galaxies (—19^ Mb ^ — 15 or 
50^ Vln/kms^^ J 120) show relatively low scatter in their SF ac- 
tivity, implying factors 2-3 fluctuations in their SFRs; above Vm ~ 
120 km/s the sequence turns off toward lower levels of SSFRs 
and larger bulge-to-disc ratios. These results are for nearby galax- 
ies, where selection effects are minimal, and imply that the SS- 
FRs of disc galaxies with Af, > 5 10* M0 follow a relatively tight 
sequence, without strong fluctuations. For galaxies smaller than 
Vra ~ 50 km/s (dwarfs) the situation seems different. The results 
by Lee et al. (2007) show that a significant fraction of such galaxies 
are undergoing strong episodic SF fluctuations due to the large scat- 
ter in their SSFRs. Another observational study of nearby galaxies 
by James et al. (2008), concluded also that there is little evidence in 
their sample of predominantly isolated field galaxies of significant 
SF through brief but intense star-burst phases. Therefore, it seems 
that the tight sequence found for normal star-forming (disc) galax- 
ies in the SSFR-Af, plane in large surveys as SDSS (Brinchman 
et al. 2004; Salim et al. 2007; Schiminovich et al. 2007) is intrin- 
sic and due to a high degree of temporal self-regulated SF within 
individual galaxies. This sequence (called the 'main sequence' in 
Noeske et al. 2007a) seems to persist back to redshifts 2 ~ 1 as 
discussed above. 

6.4 Outlook 

If the SSFR-DS phenomenon is definitively confirmed by obser- 
vations, then, as we have argued, it poses a serious difficulty for 
current disc galaxy evolution models in the context of the hierar- 
chical ACDM scenario. This difficulty, at one or another level, is 
present also in other galaxy evolution approaches, e.g., the SAMs 
(Somerville et al. 2008; Fontanot et al. 2009; Lo Faro et al. 2009). 
These models show, more from a statistical point of view than at 
the galaxy individual evolution level, that the population of rela- 
tively small galaxies (Af, ^ 10^ — lO'^'^'^ M0) is mostly assem- 
bled at higher redshifts becoming older, redder, and with (much) 
lower SSFRs at later epochs than the observed galaxies in the same 
mass range. 

Our models clearly show that less massive galaxies assemble 
their stars early because their dark halos assemble early (earlier 
on average than more massive ones), and that SN-driven outflows, 
which are more efficient as less massive is the galaxy, contribute to 
quench later SF. Gas re-accretion helps to increase moderately the 
SSFR, but not enough as to agree with observations for low-mass 
galaxies. The SF-feedback efficiencies used here to produce galac- 
tic outflows able to recover the required A/h-Af, relation are too 
large according to hydrodynamical simulations (c.f. Oppenheimer 
& Dave 2008; Dubois & Teyssier 2008). Besides, it should be taken 
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into account that a significant fraction of the SF-feedbaclc energy 
is actually dissipated into the disc ISM turbulence. Concerning the 
SF physics in our models, it should be said that the SF is assumed 
stationary (self-regulated) and related only to disc internal pro- 
cesses (isolated galaxy). Most of the model predictions (dynam- 
ics, structure, gas fractions, etc.) describe well present-day normal 
disk galaxies, but we stress that we are not modeling, for example, 
interaction-induced SF and/or SF in a bursting (non-stationary) 
regime. 

The delay of SF activity apparent from the SSFR-DS trends, 
could be produced by external effects; for example large-scale gas 
pre-heating or by the introduction of new physical ingredients in 
the intra-halo medium hydrodynamics and gas cooling process. Fi- 
nally, it could be that the SSFR-DS problem, more than a physical 
phenomenon related to the cosmic gas accretion, is a manifestation 
of some local process regarding star formation in environments that 
change with the galaxy mass and cosmic time and/or of a varying 
stellar initial mass function. 
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